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Abstract. . The dynamics of the exciton-photon entanglement, in a 
semiconductor microcavity is analyzed. Finding a closed analytical expression for 
the time evolution of the concurrence. Using as model two coupled, quantum 
oscillators with detuning between them. Supposing that the system is at a 
temperature closer to zero. And that the quantum state of the system, remains in 
the one excitation sector of the total Hilbert space. The former theoretical setup, 
describes the excitons coupled to the photons in the microcavity accurately, for 
a system with ultra low density of excitons. We remark in first place, that our 
closed expression for the concurrence dynamics. Enable us to establish a good 
analytical criterion, to determine the coupling dynamical regime, in which the 
system is. Even more, we show that our criterion, is in good agreement with 
typically accepted, experimental signatures for strong, and wake coupling regimes. 
Finally we present theoretical evidence, suggesting that the entanglement between 
the excitons and the photons. Can be measured in experiments, via the cross 
correlation contribution, to the time resolved spectra of the system. 
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1. Introduction 

Quantum information science [1], study the information processing tasks that can 
be accomplished using quantum mechanical systems. And the entanglement, is a 
uniquely quantum mechanical resource, that play a key role, in roughly all of the most 
interesting applications of quantum information science. One of the quantum system 
configurations, where is best know how to quantify the amount of entanglement, 
between two component subsystems, is a two qubits system. Where the amount of this 
physical resource, is measured by the Wotters concurrence [2], However to determine 
the concurrence experimentally. Total or partial access, to the quantum state of 
the system is quintessential. Recently various proposals had made, to achieve this 
challenge, in a two qubits system. Quantum tomography [3j is one of the most popular 
approaches, because was the first one, successfully implemented for light polarization 
degrees of freedom [4]. But have some potentially disadvantaging features, one of them 
is that have many copies of the state of the system at hand is mandatory. In order to 
improve the reliability of the results. Other disadvantage is that the information of 
system is retrieved, through local projective measurements in a quantum states basis. 
That for the most general situations could be ambiguously specified [5]. Therefore, 
great effort had made in order of improve the quantum tomography technique. Among 
them, the Ancilla- Assisted Quantum Process Tomography 6]. That proposal, reduces 
dramatically, the mandatory number of copies of the state of the system, that must 
be available, to obtain reliable results. At expenses of the coupling to a third party 
called ancilla. The direct characterization of quantum dynamics 7 . That relies 
on quantum error correction algorithms, and still having as a very hard feature. 
The unavoidable need, of have at hand a big number of copies, of the state of the 
system. In order to produce reliable information about it. All the former approaches 
has been designed to, and tested in, full optical quantum systems. Despite of that, 
some proposals were making for quantum states in other kind of architectures too, by 
example for solid state systems [8] , nonlinear optical systems j9] , and for cavity QED 
systems |10j . However nowadays, one of the most pragmatic and promising goals 
of the quantum information field, is build up interfaces between high speed flying 
qubits, as the photons used for quantum channel communications, and a stationary 
solid state qubits system. Performing quantum information processing tasks [11] . 
Material excitations in solid state systems, as the exciton for example, are promising 
candidates for solid states qubits. Because, not only offer very good coupling with the 
light field in a semiconductor microcavity -reason for which we call this component 
in our model, the matter emitter-. But also, because the exciton constitute a solid 
state system, that could provide on chip integration with other quantum computation 
hardware. The principal aim of this paper, is study the dynamics of the matter- 
field entanglement in a semiconductor microcavity at a very low temperature. By 
means the use of a "Simplified" model, of two coupled and detuned oscillators. That 
we use to provide insight to, an experimentally observable, and single copy based, 
measurement of the matter-field entanglement. Via the cross correlation spectra of 
the system. With this goal in perspective, we organize the paper as follow. In the 
second section we describe the theoretical model, its principal features, limitations 
and sketch, how the evolution in time of the state of the system, and the emission 
spectra was obtained. In the third section we present the principal results of the 
paper. Deriving a closed expression for concurrence, and showing why the cross 
correlation function, carries information about the exciton-light entanglement in the 
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microcavity. That other spectrals measurements like the photoluminescence spectra, 
cannot provide. In the fourth section we explain some features, and implications of 
our result for the concurrence. Deriving from it, in first place an analytical criterion 
based on the experimentally controllable parameters of the system, to be in the strong, 
or weak coupling regime of the dynamics. And using the photoluminescence, we show 
agreement of our criterion with often accepted signatures of the strong coupling regime. 
In second place we show the times of maximum and minimum concurrence as function 
of the experimentally controllable parameters, obtained from our closed expression for 
the concurrence. In the fifth section some concluding remarks was done. 

2. Theoretical model and dynamics of the system 

With the aim of generalize the description, to other kind of bosonic systems, we will 
refer by matter emitter, to the excitations in the material. Excitons in the ultralow 
density limit are a particular case of which. We model the system photons-excitons 
introducing a harmonic oscillator of frequency described by bosonic operators e {&) 
for the destruction (creation) of an excitation in the matter emitter, and coupling 
the former oscillator to a second one of frequency uj described by bosonic operators a 
{a)) for the destruction (creation) of a photon in the microcavity. Our interest here 
is study dynamical processes in which the energy is approximately conserved, thus 
the dynamics of the coupling between the two oscillators was taken in such way that 
rapidly rotating terms banish in the so called rotating wave approximation, we take 
the emitter-photons interaction proportional to a coupling parameter g that measure 
the rate of energy exchange between the oscillators, with the former considerations, 
the Hamiltonian of the system can be written as 

H = hfte^i + Hutfa + hg(a) e + ae f ), (1) 

where the frequency difference A = f2 — uj is often called the detuning parameter. 
With the aim of include effects of finite linewidth for the excitations in the system, 
we introduce the parameter k, that is inversely proportional to the lifetime of an 
excitation in the matter emitter. And a parameter 7, that is inversely proportional 
to the lifetime of an excitation of the field in the microcavity. Therefore, we describe 
the dynamics of the state of the system, using a phenomenological Liouvillian in the 
Lindblad form, at zero temperature and in the Markov approximation. Where each 
oscillator is initially in a coherent state 

^- = — — [H, p] + r ){2apa) — a) ap — pa) a) + n{2epe) — e^ep— pe i e).(2) 

This formulation of the dynamics have desirable properties well established in the 
literature as hermiticity, trace conservation and have semipositive eigenvalues for the 
state of the system [17] . To find the system's dynamics in ^ wc use the method 
proposed in [J5] postulating the Anzats 

p(t) = 2?(/3(t), a{t))p{t)V\p{t), a(t)), (3) 

where T>(/3(t), a(t)) is an field-matter displacement operator acting in the two modes 
with complex coherent amplitude a(t) and /3(t) for the field and the matter emitter 
respectively, and p{t) is a state of the full system in the interaction picture defined 
by the field-matter displacement operator which time derivative was reported in [15] : 
the Anzats in the equation ^ enable us to decouple the evolution of the coherent 
amplitude, from the evolution of the state of the system, finding the following system 
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of differentials equations for the evolution of the state of the system, the respective 



equations for the coherent amplitude is reported in equation ( 55 1 in the appendix 2 



m= E piling) {nh\ 



(0) 



— = — -\H, p] + r y(2apa) — a) ap — pa) a) + n(2epe) — e^ep — pe^eM) 
dt n 

Quantum effects like entanglement and purity are completely encoded in the state 
of the system p. We have choosen as initial state of the field in the microcavity, a 
coherent state p2(t) — tri/5(t) = \ato) (ao\ choosing the amplitude in such way that 
| ao |< 1 that in the case of incoherent pumping is not present is not a demanded 
restriction. In the other hand, the matter emitter is initially in a pure state that is 
linear combination of the states with zero and one excitation it is 

V (0, a ) (cos(0) |0) + sin(0) |1)) (cos(0) (0| + sin(6>) (1|) ® |0) (0| P f (0, a ) , (5) 

m 

if the initial state of the matter emitter is the excitations vacuum, the system stair in 
a state that is pure, separable and equal to p(t) = \/3(t)} {P(t)\ <g> \a(t)) (a{t)\ . Writing 
the state of the system in the base of zero and one excitation in each subsystem, it 
can be ordered as the density matrix for two qubits 

/ p88(*) pW ») o\ 

p8?(<) P° Q m Pint) o 

P i°o(t) P%(t) P\° (t) o 
\ j 

Where we truncate the Hilbert space to one excitation sector. It is we assume that 
elements of the density matrix, with more than one excitation are zero. The former 
assumption remains valid for all times in the evolution, for initial conditios with 
I Q:o lj | A) |<C 1, and the total number of excitations in the excitons-photons 
initial state equal to one almost. For this set of initial conditions, the system state 
dynamics occurs fully in the one excitation sector of the Hilber space of the system. 
Because the sectors of the density matrix with more than one excitation does not 
populate significantly. With the aim of include finite linewidth in the time evolution 
of the state of the system. We introduce the Liouville super-operator [17] 

£(•) = -I [H, •] + 7 (2a • a) - {a f a, .}) + n(2e • e f - {gtg .}). (7) 
n 

In such way that we must confront the differential equation p{t) = C(p(t)). To 
solve it, we appeal to the Laplace transform properties. Thus taking the Laplace 
transform to both sides of the differential equation -here represented as p(s)- we have 

sp(s)+C(p(s)) = p(0). That live in our hands, the simpler problem of solve a system of 

linear equations (si + £)p(s) = p(0) where the state of the system has been reordered 
by his rows as a larger vector and the Liouville super-operator as the matrix 

/ -2k -1-1 \ 

k + iSl -ig \ 



h i(u - (1) 





(8) 



ig 7 + n — i(u> - U) -ig I 

\00 00 ig -ig 2-f / 

the subsidiary condition Det((sl + £)) ^ for the existence of the solution is fulfilled, 
thus the evo lutio n of the state of the system can be find as the inverse Laplace 

transform of p(s) = (si + £) _1 /5(oj. The former is the simpler particular situation 
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of the case in which the systems is affected by coherent pumping of photons, in this 
situation we take the Hamiltonian of the system as 

H v = Mle'e + hcua^a + hg{rfe + aeJ) + hF (ae^o 1 + aV^') , (9) 

where the parameter F represents the photons pumping intensity. Switch into the 
interaction picture defined by the auxiliary Hamiltonian 'Ho — ^ D (e^e + a^a) enable 
us to write the Hamiltonian of the system for this case without time dependence 

H v = HSeJe + H(5 + A) a) a + Kg{a)e + ae 1 ) + HF (a + a 1 ) , (10) 

where a new parameter for the detuning respect to the pumping frequency S = u) — lu d 
has been introduced, and the matter-field detuning parameter A have the usual 
meaning. With the aim of include finite linewidth effects, we again describe the time 
evolution of the state of the system by means the following zero temperature Liouville 
super-operator in the Markov approximation 

l v (») = - l -[H T ,»]+ 1 {2a»a) - {a f a, •}) + n{2e • e f - {e f e,.}), (11) 

thus again the use of the Anzats in the equation |3]) enable us decouple the time 
evolution for the coherent amplitude and the time evolution of the state of the system 
with the pumping affecting just the coherent amplitudes, in this case the time evolution 
of the state is 

j t P = l v ^(p) (12) 



the coherent amplitude evolves according the equation ( 58 ) in the appendices, that 
must be solved using the initial conditions in the equation ([5j) taking again a(0) = oto 
and /3(0) = 0. The Liouvillian that we find for evolution of the system can be written 
as 

£(•) = •] + 7(2a • a f ~ {a f a, •}) + *(2e . gt _ {^e, .}), (13) 

where an effective Hamiltonian has been defined in the following way 

H V7l = Me f e + h (S + A) a) a + hg{a)e + ae f ). (14) 
Finally we remark that, the replacement rules Q, — > 5 and u) — » 8 + A relate the 



equation Q to the equation (12 1 which solution has been found before in the case 
without coherent pumping. 

2.1. First order coherence functions 

We define as a set of first order coherence functions, a set of relevant two times 
operators of the form 

G^(t,r) = (dHt)6(t + r)), (15) 

where O represents operators acting on the first or second harmonic oscillators, in our 
particular case, result of relevance the following set of coherence functions 

GW(t,T) = {tf{t)a(t + T)) (16) 

G«(t,r) = <et(t)e(i + r)) (17) 

G^ i (t,r) = (a\t)e(t + r)) (18) 

G£Ut,T) = {ei(t)a(t + T)) (19) 
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where G^(t,r) measure the contribution to the photoluminescence spectra, that is 
result of stimulated emission processes in the matter emitter, (t, r) measure the 
contribution to the photoluminescence spectra, that is result of spontaneous emission 
processes in the matter emitter, and finally (t, r) measure the cross correlations 
spectra, originated by the energy feedback between the filed in the microcavity, and 
the matter emitter. To obtain this coherence functions we appeal to the quantum 
regression theorem, the theorem states that if the evolution of the single operator 
{Oi{r)) is determined by the following equation 

^(Oi(r))=J2C j , i (d j (T)) (20) 

3 

where Cjj are a set of fixed coefficients determined by the dynamics of the system, 
thus the same set of coefficients enable us to determine the dynamics of the coherence 
functions 



d (d k (i?6i(t + r))=J2 CjAOkirfdjit + r)>. 



dT 



(21) 



In the present situation the coefficients Cjj are determined in the Heisenberg picture 
by the Liouville super-operator 



£(•) = ^ [H 7 •] + 7 (2a • a T - {a 1 a, •}) + n(2e • e f - {e T e, •}), 



(22) 



thus representing the former Liouvillian as a super-operator in the two qubits basis, 
truncated to the one excitation sector of the Hilbert space, by the matrix £, this 
matrix representation of the Liouvillian can be written in the following way 















-2k 











-27 





k + iQ 


-ig 























-ig 


7 + iui 
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ig 




















2k 


-ig 





ig 

















-ig 


7 + k + i(u> — fl) 








ig 











ig 








7 — iui 




















ig 








7 + K — i(u) — fl) 


-ig 

















ig 





-ig 


2 7 



The operators on which £ can act are represented by column vectors, as an example 



(0,0, 1,0, 0,0, 0,0,0)*, 



— 



(0,0,0,0,0,0, 1,0, 0)t, e -s- e' = 
(0,1,0,0,0,0,0,0,0)*, and e+ -> ~t* = (0,0,0,1,0,0,0,0,0)+. The former super- 
operator determines the time evolution of the relevant operators of the system, 
enabling us to obtain the coefficients involved in the quantum regression theorem 

±{a) = {it ■ £ • -*) (a(r)> + (~t • £ • it) (e(r)) = ( 7 + iu){a{r)) - ig(e{r)) (24) 
d 

oh 



(a) = {it ■ £ • It) (a(r)> + (Y • £ • ft) (e(r)) = (7 + iu){a{r)) - ig(e{r)) 
(e) = {?•£■-?) (e(r)> + (if ■ £ • ~t) (e(r)> = (« + ifi)<e(r)> - i 5 <a(T)), 



(25) 



in such way that we find straightforwardly the following set of differential equations 
for the first order coherence functions 

/ GW(t,r) \ ( 1 + 





GW{t,r) 
G^{t,r) 

\G^{t,rr J 



V 



-w 







n + ifl 


-ig 



K 



-ig 

1- in 






-ig 


7 + iw J 



.(23) 



/ GW(t,r) 
V GW(t,r) / 



.(26) 
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Defining the coefficients matrix as 21, we appeal to the Laplace Transform method 
to look for the only one solution, in such way that the coherence functions will be 
determined by the inverse Laplace transform of the matrix (21 + si) -1 acting on the 
set of initial conditions 

/ 0) \ 

G^(t,0) _ \< -\> » \')f _ /•'MM/ _ /'ni'M 

C?W (t,0) _ ~ /rt&uw ~ »oiu\ K ZI ) 

where in the intermediate step, the properties of the dynamical semigroup has been 
used [T7] . Defining the constant G = W— 4g 2 + (7 — k + i(uj — J7)) 2 we find the 
following solutions, for the set of coherence functions 







f 


(aHt)a(t)) \ 




/ (ata(t)) \ 




( Pirn \ 








(gt (t)g(t)> 




(etg(t)) 












(fit(t)e(t)) 




(ate(i)) 




p%(t) 


J 






(et(t)o(t)) y 








V pl1(t) J 



.(1) 



(t, t) = 



1 

e 2 


T(7+K 


+*( 




(p}g(«)Gcosh( 




+ Sinh 


' V 


(Pl8(*)(-7 + * " 


- iw + 


£2) + 2 


P?o(*)9)) 


1 

e 2 




+i| 




'pg}(t)Gcosh 




G*G 
+ s i„h 




(2i P j}(t) S + P g; 


(t)(T - 


« + i(« 


j - O)))) 


1 

e 2 


t( 7 +« 


+*( 




(p?J(«)G cosh ( 




G*G 
+ sinh 




(2ipJg(t)B + P?o 


(t)(7 - 




" - ")») 


1 

e 2 


t{i + k 


+*( 








G*G 
+ sinh 




(poi(')(-7 + « - 


- iu + 


Q) + 2 


PoiWs)) 



(31) 



Cr 2 ■ G*G 

According to the Wiencr-Khintchine theorem, the contribution to the spectra of any 
of the coherence functions can be obtained as 



h(t,uJ f 



1 



-me 



e lu F T Gf\t,T)d,T 



(32) 



with / taking values {C, X, Cr\, Cr2}, in our over-simplified model of coupled 
oscillators, the former integral can be evaluated analytically for each one of the 
system's coherence functions. 
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3. Results 



A remarkable benefit of the factorization of the state of the system proposed in the 
equation Is that the quantum correlations are completely encoded in the displaced 
"backward" state of the system p(t'). Here our aim is study the time evolution of the 
exciton-photon entanglement as measured by the Wooters Concurrence [2J. That for 
a displaced state of the system ordered as in the equation Q. Can be calculated 
in the following way. In first place, we define the spin-flipped state corresponding 
to p, psf = {py ® °~y)p*(°~y ® °~y) where o~y is a paulimatrix, then we calculated 
the eigenvalues of the matrix R = \J \fppsF\fp-, thus the concurrence is defined as 
C(p) = max{0, Ai — A2 — A3 — A4}, where the Aj are the eigenvalues in decresing order, 
in terms of our displaced matrix elements, the concurrence reduces to 



C(t) 



C(t) 



z % /2 g sin 2 (0)sin(^) £ 



-(7+n)i 



Fl 



(33) 
(34) 



where we has defined the following set of auxiliary constants that are functions of the 
parameters of the system 



Fl = - 7 2 



A 2 - k 2 + F2 + 



(35) 
(36) 



F2 = y / 16g i - 8g 2 (7 - A - k) (7 + A - k) + (7 - k) 2 + A 2 . 

We now will make some remarks about the units we using in all the graphics in the 
paper. The energy measured in mili-electronvolts and the time measured in units of 
-, taking 5 = 1 meV and % = 1. Thus one unit in time corresponds approximately 



to t — 6.58 x 10 fs. The expression (34) shows that the concurrence in the strong 



coupling regime, is an oscillatory function of time. That is slightly damped by effects 
of the decoherence as the times goes to infinity presenting regularly spaced in time 
collapses and revivals. We present a plot of the expression for the concurrence (34) in 
the figure |l]). 
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Time(x) 

Figure 1. Concurrence as function of time and detuning in the strong coupling 
regime ( |34[ |, with g = 1 meV, 7 = 0.1 meV and 6 = 7r/2; we choose for this 
graphic the parameter k = 0.01 meV that fulfill the conditions established in the 
equation Bgt. 
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Our main goal in this paper, is to show that an observable associated with Wotters 
Concurrence can be found in cross correlation experiments as it has been proposed 
previously for a one-atom Laser in |18j . where the cross correlation function was 
being measured, multiplying in a correlator, the heterodyne photocurrent produced by 
the photons in the cavity, by the heterodyne photocurrent produced by the photons 
generated by spontaneous emission in the matter emitter. For this purpose, now 
we will study the time resolved behavior of the quantum cross correlation function 
spectra, showing that near to a specific frequency that we will call a/ -* the matter- 
field entanglement can be measured directly from the time resolved cross correlation 
function spectra. We calculate the spectra related to the first order cross correlation 
function G' 1 ' (t, r) via the real part of the Fourier transform of this function as it has 



been defined in the equation (32 ), our interest in the quantum cross correlation is due 



to the fact that this function depends explicitly on the coherence that measure the 



entanglement, as it can be seen in the (28) and was proved in the equation (33). The 
result for the spectra calculation is 

2y F ' ir \—ig 2 + (2uj f -lu -n + i(~f + k)) 2 + (7 - k -iA) 2 J v ; 

where the information about the coherence of the d ensi ty matrix Po?(t), that 



determines the matter-fie ld c oncurrence of the equation (34) is present. Even more, 
comparing the equation ( |37[) w ith the expression for the light contribution to the 
photoluminescence spectra (|46[), just in the case of the quantum cross correlation 
function, the coefficient that multiply the coherence Poi(i) ^ s a function of ui F that 
has a local extremal value for a real valued observational frequency. Due to this 
fact, in experiments just the quantum cross correlation function enable us to extract 
information about the coherence that determines the matter-field entanglement. With 
the aim to obtain an explicit expression of the contribution to the cross correlation 



spectra of iev 2 (t, w F ), we rationalize the equation (37) that enable us to write 
I C r 2 {t: UJ F) = A P ((jJ F )pQl(t) + Ac(u F )Poi(t), where wenave defined the amplitude 
of the population and of the coherence in the following way 

A ( s ff 2 K - O) + (a, - ay) ((ay - ») 2 + k 2 ) 

C(UJf) 7r(^ + 2 S »(7K-( W ,-a;)(a;,-n)) + ((a; 1P -w) a +r ! )(( w ,_n)a + «a)) W 

a f \ = gj-ajpij + k) + + kuj) 

P{ F> n(g* + 2 5 2( 7K _ ( Wf _ w )( Wf - n)) + ((ay - o>) 2 + 7 2 ) (K - «) 2 + « 2 )) ' 1 ' 
Thus, the simpler guess for a criterion to determine the frequency where the coherence 
pio(t) could be measured from the quantum cross correlation spectra in an experiment 
as the proposed in |f 8) . is obtained maximizing the numerator of Ac{uj f ) with respect 
to uj f , the contribution of this amplitude is maximum at 

ujW = - (- V3g 2 - 3k 2 + d 1 - 2ojQ, + fi 2 + w + 2ft) , (40) 

that is always a real frequency, and therefor constitutes a good observational frequency. 
To illustrate better how the criterion works, we plot the absolute value of the two 
amplitude as function of lo f in, and out of resonance in figure ([2]), showing with 
an arrow the zone where the observational frequency uj^ is located. With the aim 
of study how the criterion works for a broad region of the parameters space, we 
have calculated the concurrence in terms of the cross coherence function contribution 

to the quantum cross correlation function spectra as 



Ac(4 0) ) 



we present this 



result in a plot in the figure ([3]), we also have calculated the error in the concurrence 
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CO p (meV) <B F (meV) 



Figure 2. Absolute value of the amplitude Aq and Ap at resonance (A = 0) in 
A and with detuning (A = 5 meV) in B. Using the parameter values g = 1 meV, 
9 = 7r/2, u) = 1 eV, re = 0.01 meV, and 7 = 0.1 meV. In the zone signaled by 
the arrows, the contribution of the coherence is always grater than the population 
contribution. 



determination by this approach as 



C(t)- 



2/cv 2 (i,c/ >) 



a plot in the figure ( 3 1 show 



this result for our oversimplified model, the guess of measure the concurrence airectly 
from the cross correlation function spectra works relatively well in a broad region of 
the relevant parameters space. It is remarkable that the result in figure ^ shows 
that the optimal behavior of the present proposal for the measure of the concurrence, 
works very well for detuning beyond 2.5 meV. It is important to remark that in 
the same spirit that we determine the criterion to measure the concurrence from 
the crossed coherence function, the maximization procedure can be done to measure 
other observable quantities from the contributions to the quantum cross correlation, 
and photoluminescence spectra. 



4. Discussion 



Based in the expression for the concurrence (34) we identify the dynamical coupling 



regimes, analyzing the four values of k that cancels the constant Fl, it is 



k= i (2 7 ± ^2^1 



2A 2 ± \/8A 2 + 1 j , 



(41) 



the former result remain valid after the roles of the parameters 7 and k are exchanged. 
The system is in critical damping regime when Fl — > 0, in this regime the concurrence 
takes the value C(t) = 2 \gt sin 2 (9)e^^ +K " >t \ . In order that the system be in the 
strong coupling dynamical regime, defined as the regime of the dynamics in which the 
eigenstates of the system are entangled matter-light states. -Fl must be a real number, 
it is fulfilled in any of the following four possible configurations of parameters 

A < -1 A k < I ( 2 7 + vW 2A 2 + y/8A 2 + 1 + 1 ) (42) 



-1<A<0Ak< 



1 



27- 



2v/2A 2 - a/8A 2 + 1 + 1 



(43) 
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0<A<1Ak<^2 7 + V2\f2A 2 - V8A 2 + 1 + 1^ (44) 

A>1Ak<^2 7 + y/2y/ 2A 2 + V8A 2 + 1 + 1^ . (45) 

Again, the former result remain valid after the roles of the parameters 7 and re are 
exchanged. It is remarkable that, in this situation the dynamical regime does not 
depend on the matter-field coupling strength g. By the other hand, the strong coupling 
regime of the dynamics is often established through the appearing of an anticrossing 
as function of the detuning in the photoluminescence spectra. With this in view, 
we will study the stimulated emission process contribution to the time integrated 
photoluminescence spectra that is defined in terms of the first order correlation 
function G^'(t,T). Introducing the constant g = ^] — Ag 2 + (7 — iA — re) 2 , we find 
the following result for this contribution to the photoluminescence spectra 

- M- ^ )g ?f ( t ( "rv:»3 ■ (46 > 

7r |_ 9 + ( 2w f + H7 + K + H w + "))) J 
To obtain the contribution to the time integrated photoluminescence spectra, the 
former result must be integrated in time as 

/cK)= lim / I c (t,uJ F )dt. (47) 

We present a plot of Ic{u F ) for some values in the detuning at the figure $5h. 
The typically accepted, signature that the system is in strong coupling regime, is 




Time(x) 

Figure 3. Concurrence as function of the detuning and the time, measured 
by the cross-coherence function in the strong coupling regime. With parameters 
a = 0.01, = it/2, 9 = 1 meV, u = 1 eV, K = 0.01 meV, 7 = 0.1 meV and 
A = 5 meV. 
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Time(x) 

Figure 4. Error in the concurrence as measured by the cross-coherence function 
in the strong coupling regime. With parameters ao = 0.01, 9 = n/2, 3 = 1 meV, 
uj = 1 eV, K = 0.01 meV, 7 = 0.1 meV and A = 5 meV. 



comprised in the apparition of an anti-crossing in the photoluminescence spectra, 
near to resonance. It can be seen in the figure ([5]), where each branch in the spectra 
corresponds to a polaritonic excitation in the system. Where our criterion for the 



strong coupling (42 1 has been used, in order to show consistency with the typically 
accepted criterion e.g. [14j . Interpolating the photoluminescence spectra profiles, for 
each value of the field-matter detuning parameter, with a Lorcntzian profile. We 
obtain the linewidth for each polariton branch, the guess of Lorentzian interpolation 
is in good approximation for our model of coupled oscillators in the ultralow density of 
excitons regime. The linewidth of the polariton branches calculated within the present 
setup can be seen in the figure ©. The crossing in the linewidth at resonance that can 
be observed in the figure ([6j constitute another signature of the strong coupling regime 
of the dynamics that is typically accepted. It is important to remark here, that the 
values in the parameters that we use corresponds to microcavities with a lower quality 
factor than the used in the reference [14] . fact that can imply a strong discrepancy in 
the states linewidth when compared with the result reported there in. By means of 



an analogous reasoning, that the one that conduce us to the result ( 42 1 , in order to 
have, the system be in the weak coupling dynamical regime, defined as the dynamical 
regime in which the eigenstates of the system are very closer to the bare matter-field 
states, Fl must be pure imaginary number. This condition is accomplished when 



1>\ } 1 VlSff 4 + 8g 2 A 2 + A 2 + k, (48) 
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Figure 5. Photoluminescence spectra as function of the field-matter detuning 
parameter, the left branch corresponds to the polariton state |xo) an d the right 
branch correspond to the polariton state \%\), see the appendix 1. Using the 
parameter values «o = 0.01, 9 = 1 meV, 8 = n/2, 10 = 1 eV, k = 0.01 meV, and 
7 = 0.1 meV. 



for all the values of the detuning parameter, again the former result remain valid after 
the roles of the parameters 7 and k are exchanged but this time the result depend on 
the value of the field matter coupling constant g. In contrast with the behavior of the 
concurrence in the strong coupling regime, in the present situation, the concurrence 
is not more periodic in time, reaching a maximum value for finite time, and tends 
asymptotically to zero as time tends to infinity, we plot this result as a function of the 
detuning and the time in the figure Q . Maximizing the expression for the concurrence 
in ( 34 1 we find the following expression for the times at which the concurrence is a 



maximum 



V2 ( 2tan- 1 ( ±V ^+2,^+2^-^-^ j + ^ 



Tm± - ^ * YT L '-, (49) 

where m run over the integer numbers and the maxima are ordered as 
7]_, t +, r 2 _, t 1+ , . . . , by the other hand, the concurrence take its minima values for 
times 

r 1 

It will result enlightening from a qualitative point of view, study the behavior of the 
maximum concurrence times, and the behavior of the concurrence at these times, with 
this aim in view, we plot thess quantities in the figure pi). Neither the expression for 
the concurrence in the equation ( 34 1 , neither the maxima or minima concurrence times 
depends on £1 and u>, depends rather on its difference A. Then in the case of coherent 
photons pumping, the concurrence evolves too as in the equation (34). 
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A(meV) 



Figure 6. In squares the linewidth as function of time for polariton branch 
corresponding to the state |xo); in circles the same for the polariton branch 
corresponding the state |xi)- With parameters ceo = 0.01, 8 = n/2, g = 1 
meV, a) = 1 eV, k = 0.01 meV, 7 = 0.1 mcV and A = 5 meV. 
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2 4 6 8 10 



Time(x) 

Figure 7. Concurrence as function of time and detuning in the weak coupling 
regime, in this graphic the time scales in units of t ~ 6.58 X 10 2 fs, with g = 1 
meV, k = 0.01 meV and 8 = 7r/2; we choose for this graphic the parameter 

7 = J 8 J_ 1 V 16 3 4 + 8g 2 A 2 + A 2 + k that fulfill the condition established in the 

equation \A&) for the weak coupling regime of the dynamics. 




i 1 1 1 1 1 1 I 1 1 1 1 1 1 1 r 



012345 012345 

A(meV) A(meV) 



Figure 8. In A, the maximum concurrence times as a function of the detuning, 
where the time scales in units of r ~ 6.58 X 10 2 fs, and the behavior of the value 
of the concurrence at this times in B, using the parameter values g = 1 meV, 
6 = tt/2, u) = 1 eV, K = 0.01 meV, and 7 = 0.1 meV. 
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Figure 9. [Color online] Tomography of the density matrix for the first four 
maximum concurrence times. With the arrow signaling the position of the 
coherence Pq5(t) in modulus for each time, the parameters we used was c*o = 0.01, 
A = 5 meV, = tt/2, g = 1 meV, u) = 1 eV, K = 0.01 meV, and 7 = 0.1 meV. 
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5. Concluding remarks 

We study the matter-field entanglement in a semiconductor microcavity. Using a 
model encompassed by two driven damped and coupled oscillators, with a detuning 
A between them. Where the first one is for the photons, and the second one is for 
the excitons in the microcavity. This Model represents accurately, a semiconductor 
microcavity, for experiments where the condition of ultra low density of excitons is 
fulfilled. It is, Na 1 ^ -C 1 where N is the number of excitons, a B is the Bohr radius 
by each exciton and d is the dimensionality of the heterostructure. We report a 
closed analytical expression that quantify the exciton-photon entanglement in terms 



of 7, k and A, using the Wotters Concurrence (34). Valid for experimental setups 
with or without coherent pumping of photons. Showing that accurated information 
about the dynamics of the concurrence can be measured indirectly. Via the cross 



correlation function (37 ). That can be measured multiplying in a correlator device, 
the heterodyne photocurrents produced by the photons in the microcavity, and the 
excitons recombination. If both are measured in a time resolved way, within a narrow 
window in frequency, centered at the observational frequency oj^ ■ 
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Appendix 1: Bare system eigenstates concurrence 

Expanding the Hamiltonian |l]) in the base with zero and one excitation, the 
Hamiltonian of the system present two eigenstates with nontrivial dynamics that are 



g |oi)-§(-A + VA2 + y)|io) 

Ixo) = , — (51) 

2 5 2 + 1 A3 - Ayjg* + f 

g |oi)-K-A- VA2 + y)|io) 

\Xi) — / \oZ) 
2gi + iA^ - A^^Tf 

often this states are called dressed states or polaritonic excitations of the system and 
each one have associated an eigenenergy 

E /i = \ (2u + A ± VA 2 + 4 5 2 ) (53) 

that define the two polaritonic branches of the photoluminescence spectra of the 
system. As it can be seen, the two states are non-separable or entangled, with the same 
Wootters concurrence -see the reference [2]- for both dressed states and determines 
the maximum matter-field entanglement available within the system, and is equal to 

Cn(A) = 2^^. (54) 
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Appendix 2: Coherent amplitude time evolution 

In the case that the pumping is not present, the dynamics of the coherent amplitude 
is described by 

d ( a \ - ( + ~ig 



n 



dt\ (3 J V -*9 -Os + iJl) J V 13 J ' (55) 

Solving the evolution equations for the coherent amplitude in the equation ( |55[ ), we 
find that they evolve in time in the following way 

e -i t [ig+y+ K +i( U +U)) ao f~ cos f gt\ + r _ K _ i(J + f Qj gin 

a(t) = ^ W ^ JJ ~, (56) 

9 

2e -t t ( 1+ K+i{ u +n)) gaQ sinff) 
P{t) = ~ (57) 

where we have defined g = — 7 — iA) 2 — 4g 2 with A — Q — uj often called the 

detuning parameter. In the case of pumping is present, the coherence amplitude 
evolves according to 

a(t) \ -ig \( a(t) \ ( uF \ (gg) 



P(t) J V ~ l 9 -K-i(S + A) J \ P(t) J \ 

that must be solved using the initial conditions in the equation ([5]), taking again 
a(0) = «o and /3(0) = 0. We find the following result for the evolution of the coherent 
amplitude 

-ItC-y+StS+iA-U) /-_ 2e I t ( T + 2 l5 + ,A + K ) FBi _ 2B2COS (£t) +2l5B3Bm (|t)] 

7 w 5 < 59 > 

((-7 + »A + K ) 2 - 4s 2 ) (g 2 t(l+ + A) + k) 2 ) 

Bj = (4g 2 - ((« - 7 ) + iA) 2 ) «J + A) - in) 

B 2 = (V - C(" - T) + »A) 2 ) (g 2 c, - < F + »!)(' - <tOK(* + A) - ««)) 

B 3 = (<*0 ((T + » 5 )C(« + A) - <») - is 2 ) (7 - >A - *) + F (2g 2 + (7 - iA - «)(<(« + A) + *))) 

1 t(7 + l(2 , + A) + K)g2 / l t(7+j(25+A) + re)pAi + f ^ (f) + am 

/3(t) = ^ Z_ 

((-7 + <A + k) 2 - 4g 2 ) (g 2 + (7 + i(25 + A) + K ) 2 ) 

Aj = (4g 2 - ((« - 7) + iA) 2 ) 

A 2 = (((* - 7) + iA) 2 - 4g 2 ) 

A 3 = (iP( 7 + «(2S + A) + k) + 2c (g 2 + (7 + iS)(i(4 + A) + *))) 
g = ^4g 2 - (7 + (iA - k)) 2 . 



(60) 
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